Some new results on one-dimensional outflow dynamics. 
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Abstract. - In this paper we introduce modified version of one-dimensional outflow dynamics 
(known as a Sznajd model) which simplifies the analytical treatment. We show that simulations 
results of the original and modified rules are exactly the same for various initial conditions. We 
obtain the analytical formula for exit probability using Kirkwood approximation and we show that 
it agrees perfectly with computer simulations in case of random initial conditions. Moreover, we 
compare our results with earlier analytical calculations obtained from renormalization group and 
from general sequential probabilistic frame introduced by Galam. Using computer simulations we 
investigate the time evolution of several correlation functions to show if Kirkwood approximation 
can be justified. Surprisingly, it occurs that Kirkwood approximation gives correct results even 
for these initial conditions for which it cannot be easily justified. 



Introduction. — The outflow dynamics was intro- 
duced to describe the opinion change in the society. The 
idea was based on the fundamental social phenomenon 
called " social validation" . By now, the opinion dynamics 
was studied by many authors, starting perhaps from the 
works by Galam [1] and developed later in the Sznajd [2] 
and Majority rule [3] models. The common feature of 
these models is that the complexity of real- world opinions 
is reduced to the minimum set of two options, + or — . 

However, in this paper we do not focus on social ap- 
plications of the model (an interested reader may resort 
to reviews [4-8]). Here we deal with a more mathemati- 
cal problem, namely finding the analytical formula for the 
probability P+ (j>) of reaching consensus on opinion + as a 
function of the initial fraction p of opinion +. This quan- 
tity is commonly called exit probability [9,10]. In fact, we 
follow the method used in [9] for the Majority-rule model. 

The one dimensional outflow dynamics is defined as 
follows: if pair of neighboring spins SiSi+i — 1 the 
the two neighbors of the pair followed its direction, i.e. 
5*1-1 — > Si{=^ Si+i) and Si+2 Si{^ Si+i); in case of 
different opinions at the central pair, the two neighbor- 
ing states are unchanged. Until now several analytical 
approaches have been proposed. One of the analytical ap- 
proaches used for the outflow dynamics was based on the 
mean field idea [11]. Within the mean field approach the 



mean relaxation time (t) as a function of the initial frac- 
tion p of opinion -|- was computed, as well as distribution 
of relaxation times. The exit probability found in this ap- 
proach is the trivial step function, at odds with the known 
simulation results for ID dynamics. 

Later, Galam in [12] presented a general sequential 
probabilistic frame (GSPF), which extended a series of 
earlier opinion dynamics models. Within his frame he was 
able to find analytic formulas for the probability p(t + l) to 
find at random an agent sharing opinion -|- at time t + l as 
a function of p{t). Among several models, he considered 
the one dimensional rule, which we investigate in this pa- 
per, i.e.: if pair of neighboring spins SiSi^i = 1 then the 
two neighbors of the pair followed its direction; in case of 
different opinions at the central pair, the two neighboring 
states are unchanged. For such a rule, within his GSPF 
calculation Galam has found the following formula [12]: 

p{t+l) = p{t)^ + lpit)'[l-p{t)] 

+ 3pmi-pit)r + ^pm-pit)r (i) 

Iterating above formula the exit probability P-^ can be 
found as a step function (see Fig. 1). It is worth to notice 
that step-like function describes exit probability in the 
case of two dimensional outflow dynamics [13], but not in 
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Fig. 1: Exit probability P+ from random initial state consisting 
of p up-spins for the modified original outflow dynamics in one 
dimension for several lattice sizes L. This is seen that results 
agree perfectly with analytical formula (solid line) given by 
eq.(23) obtained from Kirkwood approximation. Renormaliza- 
tion group (RG) results obtained in [15] for growing networks 
and calculations made by Galam within his general sequential 
probabilistic frame (GSPF) agree with simulation results much 
worse. Results obtained for modified version of outflow dynam- 
ics in which only one neighbor of central pair is changed are 
exactly the same. Results are averaged over 10'' samples. 



one dimension [14]. Step- like function for exit probability 
has been found also in [15]. 

In the paper [15] real space renormalization approach 
has been proposed to calculate the probability P+{p) of 
reaching consensus on opinion -|- as a function of the initial 
fraction p of opinion +. They have found in case of two 
sites convincing others the following analytical formulas: 

• in the case of growing network (growing hierarchical 
or the Barabasi- Albert scale-free network) 



(2) 



• in the case of fixed network they have found P+ as the 
step function observed also in computer simulations 

on the square lattice. 

It is seen on Fig. 1 that RG results for growing networks 
agree much better with simulations that RG results for 
fixed network which are exactly the same as obtained using 
GSPF by Galam. 

In this paper we present analytical results obtained us- 
ing Kirkwood approximation [9] following the method used 
in for the Majority-rule model and we obtain perfect agree- 
ment with computer simulations. Moreover, we consider 
two types of non-random initial conditions and wc show 
how analytical formulas for exit probabilities change for 
such cases. 

Approximate solution in ID. — We consider linear 
chain with L sites. Each site can be in two states ±1. We 
use the following notation: 



Fig. 2: The mean relaxation times from random initial state 
consisting of p up-spins for the modified (Inn) and original 
(2nn) outflow dynamics in one dimension for several lattice 
sizes L. In the modified version at most one spin is fiipped 
in one step while in original formulation two can be flipped 
simultaneously. Therefore in case of modifled version the time 
was rescaled by factor i. This is clearly visible the mean re- 
laxation time scales with the lattice size as ~ analogously 
to the voter model [16-18]. The results presented on the plot 
are averaged over 10'' samples. 



cr G { — 1, +1}^ state of the system. 

(j{y) state of the spin at site y if the system is in state a. 

state which differs from a by flipping spin at site x. 
Therefore a^(y) = (1 - 26^y)a{y). 

We introduce here slight modifications with respect to 
original outflow rule: choose pair of neighbors and if they 
both are in the same state, then adjust one (instead of 
two) of its neighbors (chosen randomly on left or right 
with equal probability 1/2) to the common state. Because 
this way at most one spin is flipped in one step while in 
original formulation two can be flipped simultaneously, the 
time must be rescaled by factor i. We measure the time 
so that the speed of all processes remains constant when 
L —> 00, so normally one update takes time Here, 
instead, we consider also the factor i, so single update 
takes time At — Our modification eliminates some 
correlations due to simultaneous flip of spins at distance 
3. However, if we look at later stages of the evolution , 
where typically the domains are larger than 2, simultane- 
ous flips occurs very rarely. Therefore, we do not expect 
any substantial difference. Indeed, computer simulations 
confirm our expectations - only time has to be rescaled 
(see Fig. 2). 

On the other hand, the modification simplifies the ana- 
lytical treatment. Indeed, the update rule can be equiva- 
lently formulated as follows: 



Choose randomly a spin x and side s (s = 1 for 
right, s = — 1 for left.) The updated state is 

(t(x; t + At) = a{x + s; t) if a{x + s;t) = a{x + 
2s; t), otherwise a{x; t + At) = a{x; t). 
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Within such a formulation the probability that in one 
update the flip cr — > cr^ occurs: 

+ a{x - 1)(t{x - 2) - 

- a{x){a{x + 2) + a{x + l)+a{x-l) + 



+ a{x-2)]+2 



(3) 



These flip probabilities are then inserted into master 
equation: 



P{a; t + M)^Y^ W{a' -> a)P{a'- 1) 



(4) 



Now, we make the limit L — » oo which also implies the 
continuous time limit, as At — » 0. We also note that most 
of the transition probabilities W{a' a) are zero, since 
only one spin flip is allowed in one step. Finally we end 
with 

^P(a; t)=Y^ [«^(^' ^ - wia ^ a-)P(a; t)] 

X 

(5) 

where the transition rates are trivially related to transition 
probabilities (3), w(cr^ ^ a) = 2NW{a'' a). (The sum 
is now over infinite set of sites.) For completeness we 
repeat the formula for transition rates: 



a{x + 2)(j{x + 1) + <j{x - l)a{x - 2) 



a{x) (a{x + 2) + a{x + 1) + a{x - 1) 



+ a{x-2) 



(6) 



It is hopeless to solve the master equation as it is. In- 
stead, we write evolution equations for some correlation 
functions derived from it. We define: 

Co{t) = {a{y))^Y.^{y)P{a;t) 

Ci{n;t) = {a{y)a{y + n)) 
C2{n,m;t) = {a{y - n)(j(y)a{y + m)) 
C3{n,m,l;t) = {a{y - n)a{y)a{y + m)a{y + m + 1)) 



(7) 



Only two equations are relevant for us. The first is: 



-Co{t) = -C2il,l;t) + Co{t) 



(8) 



and the second: 
d 



dt 



Ci{l;t) = -Csil, 1,1; t) - Ci(l;t) + Ci{3;t) + 1 (9) 



These two become closed set of equations, if we apply the 

approximations described in the next section. Before go- 
ing to it, it is perhaps instructive to show the intermediate 



results which lead to equations (8), (9), and analogically 
to others, for more complicated c;orrclation functions. 

So, for example, for the lowest correlation function - the 
average of one spin - we have 



dt 



{a{y)) = -2{w{a ^ a^)a(t/)> 



(10) 



and for the next one in the level of complexity 

^(^(y)^(y + 1)) = -2{wia ^ <jyMy)<j{y + 1)) 

-2{wia ^ ay+^)a{y)a{y + 1)) . (11) 

The pattern is transparent. When computing the correla- 
tion function of spins at sites xi, X2, X3, . . . , on the RHS 
we have sum of terms, in which we average the product of 
spins at sites xi, X2, xs, . . .with transition rate (which is 
constructed from the spin configuration according to (6)) 

for flip at positions xi, X2, X3, As a formula, this 

sentence means 



dt 



m^i^i)) = -2E<^('^ ^ '^"on^(^o> • (12) 



Kirkwood approximation. Now we discuss the approx- 
imations used for solving Eqs. (8) and (9). 

The first one is the usual Kirkwood approximation, or 
decoupling, which is used in various contexts and accord- 
ingly it assumes different names. For example in the clas- 
sical quantum many-body theory of electrons and phonons 
in solids, it is nothing else than the Hatree-Fock ap- 
proximation (but contrary to this theory, which may be 
improved systematically using diagrammatic techniques, 
here the systematic expansions are not developed). We use 
the name Kirkwood approximation, following the work [9] . 

In our case, the Kirkwood approximation amounts to 



in Eq. (9) and 



Cs{l,l,l;t) ^ {C,{l;t)y 



C2(l,l;t) ~Ci(l;t)Co(t) 



(13) 



(14) 



in Eq. (8). While the latter assumption (14) enables us 
to relate the equation (8) directly to (9) and therefore to 
solve it as soon as we have the solution of (9), the approx- 
imation (13) does not yet make of (9) a closed equation. 
The point is that there is also the correlation at distance 
3, the function Ci{3;t). So, we make also an additional 
approximation, which is also made in [9]. We suppose that 
Ci{n;t) only weakly depends on distance n, or else, that 
the decay of the correlations is relatively slow. If the spins 
are correlated to certain extent on distance 1 (the neigh- 
bors) , they are correlated to essentially the same extent 
also on distance 3 (next-next neighbors). This is also jus- 
tified if the domains are large enough, i. e. at later stages 
of the evolution. So, we assume 



Ci(3;t)^Ci(l;t). 



(15) 
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Fig. 3: Sample (not averaged, one simulation only) time evo- 
lution of several correlation functions given by eq. (7) for ran- 
dom initial conditions with p up-spins. Kirkwood approxima- 
tion given equations (13), (14) and assumption (15) are valid 
for later stages, although the assumption (13) agrees perfectly 
with simulation results from very beginning (left upper panel). 



In the figure 3 we present a sample (not averaged) time 
evolution of several correlation functions. Indeed, our as- 
sumptions can be justified at later stages of tlie evolution, 
although the assumption (13) agrees perfectly with sim- 
ulation results from very beginning. The second given 
by (14) agrees with simulations also quite well. Only the 
assumption (15) that the decay of the correlations is rela- 
tively slow is valid only at later stages of the evolution. 

To sum it up, the approximations (13), (14), and (15) 
say that approximately 



Co{t) 
Ci(n;t) 



m 



(16) 



where ip{t) and (j){t) satisfy the equations (dot denotes the 
time-derivative 



(1 - 0)^ 



(17) 



The solution is straightforward. We assume initial condi- 
tions (j!)(0) = mi and V'(O) = TiQ. First we solve the second 
equation from the set (17). This gives 



sinht -I- mi coshi 



cosh t + mi sinh t 
and inserting that into the first of the set (17) we have 

2mo 



(18) 



"^^^^ 1 + mi + (1 - mi) e-2* ' 
The most important result is the asymptotics 

2mo 



1 -I- mi 



(19) 



(20) 



How to interpret this finding? The average Co(t) is the 
average magnetization. In other terms, it determines the 



probability that a randomly chosen spin will have state 
+1 at time t. This probabiUty is p+{f) = (Coit) + l)/2. 
Therefore, mo = Co(0) is the initial magnetization. When 
we go to the limit t ^ 00, we know that ultimately the 
homogeneous state is reached. The asymptotic magneti- 
zation Co (00) therefore says what is the probability that 
the final state will have all spins -|-1. It is (Co(oo) + l)/2. 
So, (20) means that 



Co (00) 



2Co(0) 
1 + Ci(l;0) 



(21) 



If the initial state is completely uncorrelated, i. e. we set 
the spins at random, with the only condition that average 
magnetization is mo, we have Ci(l; 0) = mg and 



Co (00) 



2mo 
1 + ml 



(22) 



Finally, we express this result in terms of the probability 
p = (Co(0) + l)/2 to have a randomly chosen spin spin 
in state -|-1 at the beginning and the probability = 
(Co (00) -|- 1)/2 that all spins are in state -|-1 at the end. 
We have 



^^-2p2-2p+l- ^^^^ 

Computer simulations for random initial conditions, in 
which assumption Ci(l;0) = m^ can be done shows per- 
fect agreement with analytical formula (23). In the next 
section we show how results will change in case of corre- 
lated initial conditions. 

Correlated initial conditions. Here we consider 
two examples of correlated initial conditions with fraction 
p of up-spins: 

1. Ordered initial state that consists of two clusters: pL- 
length of up-spins and (1 — p)L-length of down-spins, 
for example in case of L = 10: 



p = 0.5 
p = OA 
p = 0.3 



TTTTiiiiii 
TTTiiiiiii 



(24) 



2. Correlated, completely homogeneous, initial state, i.e. 
Sn/p = 1 for n = 0, 1,2,3, for example in case of 
L = 8: 



p = 0.5 
p = 0.25 



nnnn 



(25) 
(26) 



In both cases it is easy to calculate exactly correlation 
function Ci(l;0). In the first case of ordered initial con- 
ditions we obtain: 



Ci(l;0) = l--^ 



(27) 
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Fig. 4: Exit probability P+ from ordered initial state consisting 
of p up-spiris for the outflow dynamics in one dimension for 
several lattice sizes L. Initial state consists of two clusters: pL- 
length of up-spins and (1 — p)L-length of down-spins. Results 
for original and modified rules are the same. The dependence 
between initial ratio of up-spins p and exit probability is given 
by the simplest linear function P+ = p as in the case of the 
voter model. Analytical result in this case can be obtained 
from the equation (21). In this case Ci(l; 0) = 1 — « 1 and 
we obtain Co(cx)) = co(0) = mo, i.e. P+ = p, which perfectly 
agree with simulation results. Results are averaged over 10^ 
samples. 



Thus, from equation (21): 

2Co(0) 



Co(oo) 



2mn 



1 + Ci(l;0) 1 + 1 



mo^P+= p. (28) 



Computer simulations shows that indeed for such a initial 
conditions — p (sec Fig. 4). 

As we see Kirkwood approximation surprisingly gives 
correct results also in this case. However, if we look at 
figure 5 we sec that Kirkwood approximation given equa- 
tions by (13) and (14) cannot be justified by computer 
simulations. 

We have checked also the mean relaxation time in case 
of ordered initial conditions (Fig. 6). It occurs that anal- 
ogously like for random initial conditions the mean relax- 
ation time scales with the system size as (r) ^ (sec 
Figs. 2 and 6). The same scaling has been found in the 
voter model [16-18]. However, contrary to the random 
initial conditions for which bell-shaped curve is observed, 
here the mean relaxation times is well described by simple 
parabola: 

(t) 1 

y--2P(l-f>)- (29) 

For the second correlated initial conditions, which are 
completely homogeneous we observed in computer simu- 
lations that exit probability is step like, i.e. 



antiferromagnetic state 
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for p > 0.5 
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Fig. 5: Sample (not averaged, one simulation only) time evolu- 
tion of several correlation functions given by eq. (7) for random 
ordered initial conditions with p up-spins. Initial state consists 
of two clusters: pL-length of up-spins and (1 — p)L- length of 
down-spins. Kirkwood approximation given equations(13) and 
(14) are not valid. 



In this case two-spins correlation function can be also 
calculated easily. For p = ^ < 0.5, n ^ N we obtain: 

Ci(l; 0) = p(l X (- - 2) + (-1) X 2) = p(- - 4) = 1 - 4p. 
P P 

(31) 

Thus, from equation (21): 



Co{oo) 



2Co(0) _4p-2 
1 + Ci(l;0) ~ 2-Ap 



= -1 



P+ = 0, (32) 



(30) 



which again agrees perfectly with computer simulations, 
although Kirkwood approximation cannot be easily justi- 
fied. 

Very interesting results is obtained for homogeneous ini- 
tial condition if we measure the mean relaxation time (see 
Fig. 7). Computer simulation shows that for p < 0.4 (and 
p > 0.6, respectively) the mean relaxation time (r) does 
not depend on the system size L and for p G (0.4,0.6) 
depends linearly on system size, i.e. (r) L. 

Summary. — We introduced modified version of one- 
dimensional outflow dynamics in which we choose pair of 
neighbours and if they both are in the same state, then 
adjust one (in original version both) of its neighbours (cho- 
sen randomly on left or right with equal probability 1/2) 
to the common state. We checked in computer simulations 
that accordingly to our expectations results in the case of 
modified rule arc the same as in the case of original out- 
flow dynamics, only the time must be rescaled by factor 
|. Modified version simplified the analytical treatment 
and allowed to derive the master equation. Following the 
method proposed in [9] we wrote evolution equations for 
some correlation functions and used the Kirkwood approx- 
imation. This approach allowed us to derive the analytical 
formula for final magnetization (21). In fact, just before 
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Fig. 6: The mean relaxation times for the outflow dynamics in 
one dimension for several lattice sizes L. Initial state consists 
of two clusters: pL-length of up-spins and (1 — p)L- length of 
down-spins. Results for original and modified rules axe the 
same. This is clearly visible that in case of such an ordered 
initial state the dependence between initial ratio of up-spins p 
and the mean relaxation time (t) is given by simple parabola 
not a bell-shaped curve ( like in case of the random initial 
conditions). However, stiU the mean relaxation time scales 
with the lattice size ~ L^. The results presented on the plot 
are averaged over 10* samples. 
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Fig. 7: The mean relaxation times for the outflow dynamics in 
one dimension for several lattice sizes L. Initial state consists 
of p up-spins and it is completely homogeneous, i.e. S„/p = 
1 for n = 0,1,2,3,... (e.q. for p = 1/4 the initial state is 
TiiiTiiiTiiJ.)- Results for original and modified rules are the 
same. In such a case scaling (r) ~ is not valid anymore and 
for p < 0.4(p > 0.4) the relaxation time does not depend on 
the lattice size. However, for p € (0.4, 0.6) the mean relaxation 
time scales as (t) ~ L and depends linearly on the initial ratio 
of up-spins p. The results presented on the plot are averaged 
over 10"* samples. 



finishing this paper the same result was published by Lam- 
biotte and Redner as a special case in the work [19] where 
a model interpolating the voter, Majority-rule (or Sznajd) 
and so-called vacillating voter dynamics was investigated, 
using also the Kirkwood approximation. 

In the case of random initial conditions Kirkwood ap- 
proximation can be justified looking at time evolution of 
simulated correlation functions. In this case our analytical 
results can be simplified to eq. (23) and agrees perfectly 
with simulations on contrary to earlier approaches [12, 15]. 
We have checked also how the Kirkwood approximation 
works in the case of two types of correlated initial condi- 
tions. Although in both cases the Kirkwood approxima- 
tion cannot be easily justified, surprisingly we obtained 
perfect agreement with computer simulations. 
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